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EFFECTS OF FLOOD PLAIN ENCROACHMENTS ON PEAK FLOW 

ABSTRACT 


When land development encroachments into river flood plains are permitted, 
the magnitude of the flood peak discharge will be increased due to removal of 
flood plain storage. When the flood plain encroachment is limited (for example, 
when no more than one foot of water surface rise is permitted), the study 
results indicate that the increase in the flood peak is usually relatively 
small; it is generally less than 10 percent. The arrival time of the flood 
peak can be affected to a much greater extent, and in stream systems with 
significant tributary flows, consideration must be given to the effects of 
encroachments on the coinciding of the flood peaks. Various hydrograph shapes 
and mean river slopes were used in the investigation of encroachment effects. 
Both unsteady flow simulations and steady flow profile calculations with 
hydrologic routing were used in the analyses. 


Hi 











ACKNOWLEDGMENTS 


I 

! 


The very helpful suggestions and comments provided by Arlen Feldman, 
John Peters and Bill Johnson of the Hydrologic Engineering Center are 
greatly appreciated. Typing of the report was done by Eileen Tomita. 

The report drawings were prepared by Roger Nutter, and Wayne Pearson did 
most of the computer work. Director of the Hydrologic Engineering Center 
during the study period was Bill S. Eichert. 

Funding for the study was provided by the Office of the Chief of 
Engineers, U.S. Army Corps of Engineers, Washington, D. C. 












SUMMARY 


Development of river flood plains results in removal of areas in which 
water is temporarily stored during the passage of floods. This flood plain 
storage is an important factor in the attenuation of flood peak flows, and 
also in the reduction of flood levels. As a consequence, land development 
encroachments into flood plains are usually controlled to permit only limited 
water level increases. 

It was determined in this study that when the flood plain encroachment 

is limited (for example, when the maximum water surface rise is limited to 

one foot), the increase in the magnitude of the flood peak is usually small, 
in most cases it is less than 10 percent. The arrival time of the flood peak 
can be affected to a much greater extent, however. Therefore, in stream 
systems with significant tributary flows, consideration must be given to the 

effects of encroachments on the combining of the flows, since a substantial 

increase in peak flow may result at certain locations. The analyses made in 
this study provide some general indications of the effects of encroachments 
on flood characteristics; however, the problem is such that each individual 
case must be analyzed separately if numerical results are required. Procedures 
for this are available and should be used when it is expected that an increase 
in flood peak flow would cause problems. A major purpose of this report is to 
provide guidance to help identify those hydraulic situations where significant 
increases in flood discharge peak can be anticipated. 

Encroachment effects were determined for a simple river channel and 
flood plain reach with constant slope, roughness, and width. Several flood 
plain widths were used as well as two different river slopes. Some results 
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are given for steady flow profile analyses with attenuation effects determined 
using hydrologic routing techniques. 


Unsteady flow simulation with the "Gradually Varied Unsteady Flow" computer 
program was used for the major portion of the work. Five different hydrograph 
shapes were used. It was found in this study that the effect of a given 
encroachment is less in steeper rivers than in flatter ones. Also, the effects 
of hydrograph shape are relatively independent of valley width for steeper 
slopes. The hydrograph shape is an important factor in flood wave attenuation 
for flatter river slopes, however. 

In general, the situations which produced the greatest increases in peak 
flow were those with the flattest slopes, with the largest reduction in 
storage, and with the hydrograph peak early in the flood event. If a 
combination of these factors exists, increases in peak flow can be substantial. 
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1. INTRODUCTION 


1.1 Purpose of Stud y 

The purpose of this study is to provide guidelines for the assessment of 
the effects of flood plain encroachments on flood peak flow, maximum water 
surface elevations, and flood-peak arrival times. Because each river presents 
a unique hydraulic problem due to its unique geometry, types of lateral inflow 
conditions, and other conditions, each situation must be analyzed separately. 
However, as illustrated in this report, some general guidelines can be developed 
which will indicate when potential problems can arise. If problems are indicated, 
then further study can be made to provide more definitive answers. 

Much attention has been devoted in recent years to establishing limits to 
encroachments into flood plains to reduce damages due to flooding. The National 
Flood Insurance Program (Federal Insurance Administration, 1977) has provided 
the impetus for many studies of this nature. In these investigations, the 
maximum floodway encroachments are defined on the basis of some maximum rise 
in the flood water level (usually one foot or less). However, the effect of 
the encroachments on the magnitude of flood peak flow has frequently been ingored 
in these studies, primarily because the steady-flow analysis commonly used does 
not provide a means of determining the effect of encroachments on the discharge. 

The rate of movement of the flood wave down the river can also be affected 
significantly by encroachments, as a result of the elimination of areas for 
water storage on the flood plain. However, the commonly used steady-state 
hydraulic analyses provide little information regarding wave travel times, 
other than what can be inferred from mean flow velocities, and an unsteady flow 
model may be required if flood wave movement is of interest. 










Why should one be concerned about flood peak flow values or changes in 
the timing of the flood peak? The answer to this question is: if flood 
peak flows are significantly increased by flood plain encroachments, down¬ 
stream areas may have reduced protection from floods. The flood-peak 
arrival time may also be of concern, because encroachments usually 
increase the flood-wave speed. This causes peak flows to occur earlier, thus 
giving less time to respond to flood warnings. On streams with tributaries, 
a change in the timing of the peak may be beneficial or detrimental, depending 
on whether or not the peak flows of the tributary and main-stem coincide. 

This report discusses some of the types of flood events for which 
significant increases in flood peaks may occur as a result of encroachments 
into the flood plain. It provides guidance for the engineer making flood 
plain investigations by indicating the type of situations for which additional 
effort should be devoted to computing the flow through the affected reach of 
river by flood routing techniques. This study is concerned primarily with 
situations in which the flood plain encroachments follow the criteria for 
flood insurance studies, i.e., encroachments that cause a one foot (or less) 
rise in the water surface. Encroachments causing large water surface changes 
will affect the flood hydrograph to a much greater degree and different 
procedures may have to be used for these cases. 

1.2 Study Background 

This study was initiated in response to a request from the Flood Plain 
Management Services Branch, Office of the Chief of Engineers, II. S. Army Corps 
of Engineers. An early phase of this investigation involved the evaluation of 
some typical encroachment problems. Analyses were made using computer programs 
HEC-1 and HEC-2. Flood hydrographs were determined by basin modeling with HEC-1 
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(Hydrologic Engineering Center, 1973). The flood peak flows from the HEC-1 
results were used to determine encroachment limits such that the rise in the 
water surface with the encroachments was no more than one foot. Standard 
techniques for determining encroachment limits were followed using the 
encroachment options of computer program HEC-2 (Hydrologic Engineering Center, 
1979). Several representative watersheds and flood channels were investigated 
in this work. A summary of the results of this initial work is included in 
Chapter 3. 

In a later phase of the study, the major effort was concentrated on the 
use of unsteady-flow-analysis techniques to develop relationships between 
parameters describing the flood plain and encroachments and the flood hydro¬ 
graph. A number of the earlier HEC-2 runs were duplicated using an unsteady 
open-channel-flow computer program. Simplified flood routing techniques were 
also used, and these simplified procedures also proved useful in the development 
of guidelines. The results of these analyses are described in Chapter 3. 

1.3 Literature Review 

A number of reports and papers dealing with the effects of flood plains on 
flood hydrographs were reviewed. The problem is of world-wide concern; however, 
no simple and widely applicable answer to the problem was found in the literature. 
The approach generally conceded to be the best way of evaluating attenuation 
problems is to perform an unsteady flow analysis of the particular river, using 
a computer program which solves the full unsteady flow equations, if needed. 

In most practical applications however, hydrologic or hydraulic flood routing 
techniques are used to route flood flows, and water surface elevations are 
then determined by steady open-channel flow computations using the appropriate 
flow values. Physical models are also used to determine encroachment effects. 
Although costly, they permit details of the flow patterns to be studied. 











Most of the literature reviewed dealt with the general effects of flood 
plains on flood hydrographs, and did not specifically treat flood plain 
encroachments. A comprehensive presentation of the general subject of flood 
routing in rivers is that of the Natural Environment Research Council of 
Great Britain (1975), which covers flood routing methods, theory, and 
applications. Although the material is slanted toward use with British rivers, 
it is general enough to be applicable to a wide range of river types. 

Abbott, et aL. (1971) compared the results obtained using a numerical 
model of a river and flood plain system with data from the Danube River and 
from a physical model of the Danube. They concluded that a one-dimensional 
numerical model can give good simulation of flood events if storage capacity 
as a function of stage and the effective flood plain roughness are accurately 
specified. 

Some observations of flood plain influences on flood waves for Russian 
rivers are reported by Grashevsky (1967). He also made numerical analyses of 
simple channels with flood plains. His studies indicated that flood plain 
"accessibility," i.e., the relative ability of the water to leave the river 
and enter the flood plain, was an important factor in flood wave attenuation. 

If a flood wave moves past areas in which large amounts of potential flood 
plain storage exist, but the movement of water into these areas is restricted 
by high hydraulic roughness in overbank areas or by limited openings to the 
flood plains, then the full available storage may not be utilized. Grashevsky 
also observed that in channels with slopes on the order of 0.001 the attenuation 
of the flood wave was observed to be much less than with channels with flatter 
slopes. 
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Huang and Gaynor (1977) studied the downstream effects of stream channel 
improvements on floods using a computer program which utilized kinematic 
routing. Their general conclusion was that "channel improvements" cause down¬ 
stream flood hydrographs to have higher peaks and also cause the peak flow to 
occur earlier. 

Johnson and Senter (1977) evaluated the loss of flood plain storage on 
the Ohio River for long river reaches. They simulated different amounts of 
storage loss ranging from complete removal of valley storage by use of levees 
to partial loss of valley storage with levees on only the downstream portion 
of the river (the lower 58 miles). Also examined were the effects of increasing 
or decreasing the Manning's 'n‘ values by 20 percent. They concluded that for 
the reach of the Ohio River they were studying, the flood height would not be 
substantially changed by complete loss of valley storage or by substantial 
variation in Manning's 'n'. 

Effects of proposed encroachments in the flood plain of the Connecticut 
River were studied by Dewey and Kropper (1964). They examined the consequences 
of removing major volumes of valley storage along 40 miles of the Connecticut 
River. Storage routing techniques were used in these analyses. The general 
effects noted included: 

(1) A 10 to 40 percent increase in peak discharge downstream as a 
result of the loss of 10 to 30 percent of the valley storage; 
and 

(2) An increase in water surface elevations for major flood 
discharges as a consequence of the decrease in waterway cross- 
sectional areas. 

The approximate increases in water surface elevation were: (a) 1 foot for 

a 10 percent reduction in valley storage, (b) 4 foot increase for 20 percent 
reduction, and (c) 7 feet for a 30 percent reduction. 
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A large number of hydraulic model studies involving the effects of levees 
and other types of encroachments on river water levels have been made by the 
U. S. Army Corps of Engineers, Waterways Experiment Station (WES). In these 
studies, major changes to river flow areas and flow patterns were examined. 

These frequently result in water level changes of several feet for major 
flood events. Examples of these types of studies are given in the reports 
describing the various WES investigations using the Mississippi River Basin 
Model (e.g.. Waterways Experiment Station, 1955a; 1955b; 1959; 1971). 

The effects of levees and encroachments on design flood profiles were determined 
in this very large physical model using hypothetical flood waves (for "normal 
floods" and flash floods). In addition to these unsteady flow analyses, steady- 
flow tests were also made. Stage-discharge hydrographs were measured for the 
unsteady flows, and water surface profiles for the steady flows were determined 
for two cases: before and after levees were positioned in the model. 

A comparison of various approaches is given by a study by Huntington (1974) 
which involved both a physical and a mathematical model. He found that the 
mathematical model gave results that matched the physical model quite well in 
general, although in the flood plain itself the mathematical model consistently 
overestimated the flood levels. This was considered to be due primarily to 
the ignoring of the fine-scale details of the river and flood plain geometry 
in the mathematical model as well as neglecting the head loss factors at the 
entrance to and exit from the flood plain. However, these were minor errors, 
and Huntington concluded that mathematical models of unsteady river flow can 
provide good simulation of both flood flow magnitudes and stages in the analysis 
of rivers with flood plains. 
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An additional factor that should be considered is that the actual water 
surface rise will often be less than the computed because the full proposed 
encroachments frequently do not occur. In a recent survey by Goddard (1979), 
it was found that the actual development on flood plains resulted in an 
average "flow" blockage of only 25 percent (although the blockage varied from 
zero to 100 percent). In most cases, therefore, significant flood plain 
storage is still available even though the floodway fringe is developed. 

Goddard's survey also showed that actual floodway encroachments usually 
result in a water surface increase that is significantly lower than 1.0 ft. 

For the floodways which he evaluated the mean water surface increase was about 
0.7 feet. 
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2. APPROACHES USED IN STUDY 

2.1 Flood Wave Movement in Rivers 

As a flood wave moves downstream in a river, it is perceived at a given 
point in the river as first a rise and then, after the peak has passed, as a 
lowering of the water surface. In all channels, as the magnitude of flow 
increases, the stage also increases, and water in temporary storage in the 
river will increase as well. Because some of the water entering a given reach 
of the river goes into storage, the rate of flow leaving the reach is decreased 
as the consequence of maintaining a balance between the volumes of water 
entering and leaving the reach and the volume stored in the reach at any 
given time. This decrease in flow is called the "attenuation" of the flood 
wave and occurs to some degree in all channels; Figure 2.1 shows a typical 
situation. In rivers with wide flood plains the volume of overbank storage 
can be quite large as is indicated in Figure 2.2. This storage can be an 
important factor in reducing flood flows downstream. 

As the flow decreases and the river stage falls, the water stored in the 
channel and in the overbank areas is released from storage. This lengthens 
the time base of the flood. Because the energy gradient is steeper on the 
rising side of a flood wave than on the falling side, the flow is greater on 
the rising side for a given stage. This is illustrated by the looped stage- 
discharge relationship in Figure 2.3. 

An important effect of storage on the flood hydrograph, therefore, is 
attenuation of the peak discharge. Removal of flood plain areas by fills, 
levees, buildings, road embankments, or other encroachments will reduce the 
volume of temporary flood storage, and thus there will be less attenuation 
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Flood Hydrograph Attenuation 

















I 


I VALLEY SIDESLOPES| 
^ FLOODPLAIN k. 




CHANNEL 



y 

'J- ’ . ... ( v v p 

■ "* 


/ 




/ 

TYPICAL VALLEY SEGMENT 


FLOODPLAIN 

- m 

CHANNEL 

+ ■ m 

« - 

FLOODPLAIN 


WATER SURFACE 
OVERBANK STAGE 


WATER SURFACE 
BANKFULL STAGE 




He. 


CHANNEL BOTTOM 


TYPICAL VALLEY SECTION 


FIGURE 2.2 Schematic Diagram of River with Floodplain 
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FIGURE 2.3 Typical Stage-Discharge Relationship for Rivers 
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of the flood peak. The resulting higher flows may significantly change flood 
flow-frequency relationships downstream from the encroached reach from those 
previously existing. The higher discharges may cause significant increases 
in water surface elevations. And, as mentioned above, problems related to 
changes in timing of the flood peaks may arise if combinations with tributary 
flows are significant. 

An encroachment can cause backwater effects (and hence an increase in 
channel and overbank storage) upstream from the encroachment. This can be 
as significant as the loss in storage due to the encroachment in some cases. 

2.2 Formulation of the Problem 

In studies made to establish encroachment limits (such as flood insurance 
studies), the peak flow for a given level flood, say the 1% (100-year) flood, 
is used with a steady-flow analysis to determine the floodway limits for which 
the water surface does not rise more than a specified amount (usually one foot) 
above the natural-condition profile (Goddard, 1979). In the steady flow analysis, 
the width of the waterway (and thus also its ability to carry the flood flow) is 
reduced by the amount which gives the specified water level rise. However, in 
most cases the effect of the loss of flood plain storage on the magnitude of the 
discharge is ignored in the steady flow analysis. 

For flood events resulting from intense, relatively short-duration storms, 
the flood hydrograph generally has a sharp peak. For this type of event, the 
effects of flood plain storage can be very significant. In contrast, when the 
rainfall duration is long relative to basin response time or when the flood 
results from snowmelt such that high flows persist for a long time, the rate of 
change in flow is relatively small, and steady-state assumptions usually give 
good results. 
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The engineer or hydrologist who is involved in determining flood plain 
encroachment limits usually has the following information with which to work: 
(1) either a flood hydrograph for the reach of river under study, or more 
likely, the magnitude of the flood peak, (2) a detailed description of the 
river and flood-plain geometry (such as cross sections), and (3) some measure 
of the hydraulic roughness characteristics of the river (such as Manning's 'n' 
values). The objective of the study is to define the allowable encroachment 
limits and determine water surface profiles with and without encroachments. 

When performing an encroachment study, the hydrologist should consider 
the following questions: 

(1) Will there be significant changes to the flood hydrograph due to 
the encroachment? Most importantly, will the peak flow at a given 
downstream location increase as a result of the encroachment? 

(2) How is the arrival time of the flood peak affected by the encroach¬ 
ment? 

Answers to these questions are attempted in two ways in this report: 

In the first, analytical approaches based on unsteady flow computations are 
used; in the second, information from observed flood events and results of 
hydraulic model studies are evaluated to determine the usefulness of various 
analytical approaches. 

It is assumed here that the flood hydrograph is known at some point 
upstream from the proposed encroachment. The problem therefore, is one of 
routing this flood hydrograph to some downstream point to determine not only 
the maximum flow at that point, but also the maximum water surface elevation 
associated with the flow. The problem is basically an unsteady open-channel- 
hydraulics problem. 
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There are a number of methods available for making flood routing computations 


(Henderson, 1966; Strelkoff, 1969). The basic equations for the description of 
unsteady open channel flow (the St. Venant equations) have been known since the 
latter part of the nineteenth century, but it is only since the advent of large 
digital computers that numerical solutions to the full equations of motion can be 
made. Solutions of the full St. Venant equations are not routinely used, because 
they are much more difficult and costly than other methods. Also, the available 
procedures require knowledge of both hydraulic theory and numerical analysis 
techniques, as well as access to a computer. As a consequence, most flood 
routing is done using simplified routing methods (Weinmann and Laurenson, 1979). 
In many cases, in fact, simplified methods give results that are nearly as 
accurate as those obtained from the full St. Venant equations. For example, 
for rivers with slopes greater than 0.001 some of the terms in the full equations 
can be neglected (Henderson, 1966; Natural Environment Research Council, 1975), 
and it is not necessary to use the full equations. 

2.3 Unsteady Flow Analyses 

Full Equations of Motion . The computer program which was used in this 
study to analyze the effects of encroachments on floods is titled "Gradually 
Varied Unsteady Flow Profiles, Flood Routing by the Hydraulic Method." The 
program is available from and maintained by the Hydrologic Engineering Center 
(Hydrologic Engineering Center, 1977). It is based on the "TVA Model" described 
by Garrison, et aK (1969). 

The full equations describing open channel flow are used in this computer 
program, and thus this analysis simulates the effect on the flow resulting from 
changes in flow cross-sectional area as the flood wave moves along the river. 







However, the program, because it is a "one-dimensional" model, takes into 
account only the motion of the water in the direction of the channel axis. 

It does not make provision for water movement from the channel out into the 
flood plain and back into the channel again, except in the accounting of total 
water volume. Even with these limitations however, this type of model has been 
shown to give good simulation of actual flood events in many cases. This 
approach was selected for analyses described subsequently. 

Detailed treatments of the hydraulics of rivers and flood plain interactions 
are described in the literature. Liggett (1968), for example, treated the 
flood plain as an offstream storage element that contributed only partially 
to through-flow. This type of formulation is useful when evaluating the 
storage effects of water backed-up in a tributary, for example. This can 
be handled in the Gradually Varied Unsteady Flow program by the use of "non¬ 
conveyance" areas in a cross section. 

Fread (1976) treated the flow in the main channel separately from the 
flood plain flow (both are considered one-dimensional flows). The effects 
of a meandering river in a wide flood plain can be taken into account in this 
approach since the quantity of flow as well as the flow path lengths and 
hydraulic properties are accounted for separately. 

The process of energy transfer between the main channel and the flood 
plain was considered in a flood routing approach by Radojkovic (1976). The 
effects of channel shape on resistance in flood plain channels was studied by 
Yen and Overton (1973), Wright and Carstens (1970), and Malhotra (1968). In 
these studies, the flow was treated as one-dimensional, a reasonable 
assumption for both the main channel and flood plain flows if extensive 
meandering does not occur, as indicated by the experimental study of Toebes 
and Sookey (1967). 







A simple and uniform river-channel and flood-plain geometry was used in 
this study to provide cross sections which were easy to describe. The main 
concern here is to be able to define the loss of flood plain storage due to 
encroachments in relation to other parameters. The river cross-sectional shape 
used is shown in Figure 2.4. It is symmetrical about the channel center line. 
Constant channel slopes of 0.0001 and 0.001 (0.5 ft/mi and 5 ft/mi) were 
used. A single Manning's 'n 1 was used to characterize the roughness of the 
full cross section. 

A 5.4-mile long channel reach was used as the study reach. A flood hydrograph 

1 

was introduced at its upstream end, and the flow was routed downstream. The 
downstream boundary condition was taken as normal depth plus 0.1 feet at a 
point 0.7 miles downstream from the end of the study reach (i.e., 6.1 miles 
from the upstream end). This was found to give a reasonable representation of 
downstream conditions. Runs made using the same boundary condition one mile 
further downstream gave similar results. 

The flood hydrograph at the upstream end of the reach was represented as 
triangular in form for simplicity. More critical than the precise shape of 
the hydrograph is the timing of the peak flow relative to the beginning of 
the hydrograph. To evaluate this, several different hydrograph shapes, each 
with different times to flood peak were used. The various hydrographs are 
shown in Figure 2.5. 

Simplified Flood Routing Methods . As mentioned above, in many cases, it 
is not necessary to use the full equations of motion for analyzing unsteady 
river flows. In fact, most flood routing in practice is done either by 
"hydrologic" (storage-routing) methods or by simplified "hydraulic" routing 
methods (such as the kinematic wave or diffusion equations) which are based 
on simplified forms of the equations of motion. 
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Perhaps the most widely-used of the hydrologic methods are the Muskingum 
and modified-Puls methods. In the hydraulic routing group are the kinematic 
wave method and solutions to the "convection-diffusion" equation (Natural 
Environment Research Council, 1975). The convection-diffusion equation was 
also used in this study to illustrate how various factors, such as hydrograph 
shape, affect flood wave attenuation. Additional discussion of the convection- 
diffusion equation approach is given in the Appendix. 

2.4 Search for Field Data 

A search for field data pertinent to the problems considered in this study 
was of very limited success. Flood plain management units in various 
Corps of Engineers Districts around the country were contacted. They were 
asked for information on rivers for which data might be available on floods 
which occurred both before and after encroachments (such as levees) were 
placed in the flood plain. As a result of this search several rivers were 
suggested for evaluation. However, no numerical evaluation of encroachment 
effects could be made, because sufficient data to properly describe the problem 
were not available in all cases examined. 

A major obstacle to acquiring field data that describe encroachment effects 
is that most flood plain encroachments are in areas that also experience rapid 
development (such as urbanization) in major portions of the watershed. As a 

consequence, the hydroloqic characteristics of the watershed may change greatly, 
and it is not possible to separate the encroachment effects from influences of 
other watershed changes on a flood hydrograph. 
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3. STUDY RESULTS 


3.1 General 

The natural attenuation (reduction in amplitude) of a flood wave is dependent 
on a number of factors. These include: 

(1) The area available in the overbank regions to provide off- 
stream storage. 

(2) The slope of the river. 

(3) The shape of the hydrograph, especially the curvature at 
the peak. 

(4) The travel time (or speed) of the flood wave. 

This study is concerned primarily with the effects of encroachment which produce 
a limited change in the flood elevation (encroachments which follow the guidelines 
established by the Flood Insurance Program). Therefore, the main emphasis here 
is on the evaluation of the effects of encroachments producing a one-foot rise in 
water surface elevation. Encroachments of greater magnitude which will cause 
more pronounced changes are not considered here. 

3.2 Background on Case Studies 

A hypothetical watershed and downstream channel reach (shown schematically 
in Figure 3.1) were used in the first phases of the study. This watershed and 
channel were chosen to be representative of conditions that might occur in a 
moderate-size basin (400 sq mi). An inflow hydrograph for the upper end of 
the routing channel was based on a design storm calculated using the unit 
hydrograph computational procedures available in computer program HEC-1 (Hydro- 
logic Engineering Center, 1973). Various basin slopes and channel lengths 
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were used in the computer runs to evaluate their effect on basin response 
(i.e., on hydrograph characteristics). 

Two different channel slopes and several flood plain widths were used with 
the flood plain and channel cross section shown in Figure 3.2. Both the channel 
slope and cross section were assumed to be constant along the length of the 
channel. A steady flow simulation for each individual channel width was made 
using computer program HEC-2 (Hydrologic Engineering Center, 1979). The HEC-2 
runs were made with the option which permits automatic determination of encroach¬ 
ment limits to match a target value of water surface elevation rise. The 
encroachment which caused a water surface elevation exactly one foot higher 
than under original conditions was determined for each case. The peak flow 
for the computed hydrograph at the upstream end of the reach was used as the 
discharge through the complete reach for the determination of encroachment 
limits, a commonly used procedure for these types of studies. The results 
of these calculations are presented in Table 3.1. The reduction in river top 
width for the cases analyzed ranged from 11 to 38 percent, with corresponding 
losses in valley storage of 5 to 17 percent. 

The original hydrographs were routed through five-mile and ten-mile 
channel reaches usinq storaqe-routing procedures (modified-Puls method) of 
program HEC-1. The largest increases in peak flow were about 16 percent for 
the ten-mile reach and about 9 percent for the five-mile reach. Most of the 
cases examined had much smaller peak flow increases. The situations 
which produced the greatest increases were those with the flattest slopes, 
the largest reduction in storage, and an early-peaked hydrograph. If a 
combination of these factors exists, the increases can become quite significant. 
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EFFECTS OF ENCROACHMENTS USING HYDROLOGIC ROUTING METHODS 



Increase in Stage (ft) 













3.3 Unsteady Flow Analyses 

Effects of encroachments on the magnitude and timing of the flood peak were 
determined through unsteady flow analysis using the "Gradually Varied Unsteady 
Flow Profiles" computer program. Hydrographs were routed down a 6.06 mile long 
channel, using a stage-discharge relationship at the downstream end to define 
the downstream boundary condition. The index point for hydrograph comparison 
was 5.4 miles from the upstream end of the channel. This location was judged 
to be far enough away from the downstream boundary to not be significantly 
affected by it. 

The flood hydrographs were triangular in shape for simplicity. A comparison 
between an original hydrograph as computed by the HEC-1 computer program and its 
triangular approximation is shown in Figure 3.3. This hydrograph is Hydrograph 
3 of Figure 2.5. 

Tabulation of Results . Table 3.2 gives the results of the unsteady flow 
calculations for the five hydrographs of Figure 2.5. The initial flood plain 
widths were 350, 600, 800, 1200, and 5000 ft, while the corresponding widths 
giving a 1.00 ft rise for steady flow were 292, 490, 652, 960, and 3100 ft. 

In general, the increase in maximum flood stage computed by the unsteady flow 
analysis is close to one foot, except for those cases where the flood plain 
volume is large with respect to the volume of the flood. In the latter case, 
the increase is usually less than one foot. 

Discussions of other aspects of the results are given in the following 
paragraphs. 

Effects of Hydroqraph Shape . The effect of hydrograph shape on flood 
attenuation was examined in this study using the five triangular hydrographs 
shown in Figure 2.5. The relative flood peak magnitude, Q p /Q 0 * is given in 
Figures 3.4 and 3.5 as a function of valley width W for each hydrograph 
(Qp = peak flow at point 5.4 miles from upstream end of channel; Q q = peak 
flow at upstream end). 
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Table 3.2 

RESULTS OF UNSTEADY FLOW ANALYSIS 
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The plot in Figure 3.4 shows that in narrow valleys the hydrograph shape is 
an important factor in flood wave attenuation for the flatter river slope 
(S = 0.0001). When the peak occurs early in the flood period, the attenuation 
is greater than when the peak occurs later in the flood. In the case of the 
early peak, flow from the region of the hydrograph peak supplies the volume to 
fill the flood plain storage, while in the delayed-peak case, a good portion of 
the flood plain storage has already been filled when the peak of the flood 
arrives. In wide valleys with the flatter slope (S = 0.0001) hydrograph shape 
becomes less important. 

For the river channels with the steeper slope (S = 0.001) the effects of 
hydrograph shape are relatively independent of the valley width. Hydrographs 
with early peaks experience more attenuation than those with late peaks for 
all valley widths simulated. The sample routed hydrographs shown in Figure 3.6 
for Hydrograph 1 (early peak) and Hydrograph 5 (late peak) illustrate the 
difference in the transformation of the hydrographs for the two cases for a 
valley width of 292 ft. 

Increase in Flood Peak Flow . When the water surface rise was limited to 
one foot the increase in flood peak flows was usually less than 10 percent, 
except where the flood plain was initially very wide and the encroachment 
large. In these cases the increase in peak flow can be quite large. 

Figure 3.7 shows the increase in peak flow for the routed hydrographs for 
the case where the encroachment causes a 1.00 ft rise in the water level. 

The increase is much less for a given reduction in flood plain area for the 
steeper slope (S = 0.001) than the flatter slope (S = 0.0001). To keep the 
increase in peak flow at 10 percent or less, the reduction in flood plain area 
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FIGURE 3.7 Increase in Flood Peak Flow as a Function of the 
Reduction in Flood Plain Area 






should be limited to approximately 15 percent for the flatter slopes, while much 
greater amounts of the flood plain can be occupied for steeper slopes for the 
same increase in peak flow. 

Effect of Channel Slope . The effect of channel slope on the routed flood 
hydrograph is shown in Figures 3.8 and 3.9. Figure 3.8 illustrates the effect 
of the slope on hydrograph shape for the case of a 292 ft wide channel using 
Hydrograph 1 of Figure 2.5. In the steeper channel (S = 0.001) the attenuation 
is much less, the arrival time of the peak sooner, and the flow subsides much 
more rapidly than for the case when the slope is flatter (S = 0.0001). Figure 
3.9 indicates the relationship between computed flood peak flow and valley 
width for the two slopes used in this study. The hydrograph used in this 
comparison is Hydrograph 3 (Figure 2.5). For the cases studied here, the 
hydrograph for the flatter slope was attenuated less. 

The maximum water depth on the flood plain also depends on the river slope 
as well as valley width. Figure 3.10 gives a plot of the maximum water depth 
above the local river bottom as a function of valley width for the two river 
slopes used in this study. The data for Figure 3.10 are from the same computer 
simulations as used for the data for Figure 3.9. For this hydrograph, the 
water depths on the flood plain are greater for flatter slopes and the narrower 
valleys, while for wider valleys the depths are greater for the steeper slopes. 

Another indication of the influence of the channel slope on the behavior 
of the flood wave is given in Figure 3.11, which shows the variation in water 
level in the river with time for the two slopes: S = 0.001 and S = 0.0001. 

The ordinates of the two curves are plotted relative to the original water surface 
elevation (the depth for the case of the 0.0001 slope is about 3.2 ft greater at 
the beginning of the hydrograph). For this particular case, the total rise in 
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IGURE 3.11 Effect of Channel Slope on Stage Hydroqraoh 












water level is about the same for the two slopes; however, this will not always 
be the case. The stage hydrographs for the steeper river peak earlier and 
subside more rapidly in a manner similar to the discharge hydrographs of 
Figure 3.8. For flatter slopes, water remains in storage on the flood plain 
for much greater lengths of time. 

Flood Peak Timing . The timing of the flood peak is strongly affected by 
the degree of encroachment. In general, the flood peak arrives significantly 
earlier as a result of a reduction of flood plain storage. If tributary flow 
or high local inflows are important for the stream being studied, a change in 
timing of the peak could have major consequences. If it is expected that this 
will happen, an analysis of this effect should be made by routing and combining 
hydrographs. 

The effect of valley width on the arrival time of the flood peak is shown 
in Figure 3.12 for Hydrograph 3 for two river slopes. In narrow valleys the 
slope does not appear to influence the peak timing greatly. For wider valleys, 
the flood peak arrives sooner in the steeper slope case, as would be expected. 
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4. APPLICATION OF RESULTS 


4.1 Procedures and Guidelines 

Each river and its flood plain present a unique situation, and at this point 
no specific guidelines for evaluating individual cases can be stated. Each 
encroachment problem must be considered by itself. 

However, within the range of conditions prescribed by the Flood Insurance 
Administration's guidelines (a one-foot maximum rise in water surface), it is 
possible to draw some general conclusions. As shown by the unsteady flow 
analyses discussed in the previous section, in many cases when the rise in water 
level is restricted to no more than one foot, the increase in the flood peak flow 
is small - frequently less than ten percent. When this relatively minor 
difference in flow is considered in the light of the uncertainties which 
are present in other aspects of the hydrologic and hydraulic descriptions 
of the problem (such as the flood plain geometry, roughness parameters, and 
the precise value of the flow magnitude) it seems reasonable to ignore the 
increase in peak discharge. 

The exceptions to this general rule are for the cases of rivers with very 
flat slopes and with encroachments which remove large amounts of the flood 
plain area. In these cases there can be a substantial increase in the flood 
peak magnitude. Analysis of the changed hydraulic conditions must be made in 
these situations because the downstream hydrologic conditions (the discharge- 
frequency relationship) will be changed. In most cases, standard flood routing 
techniques can be used to compute a new hydrograph for the downstream end of 
the reach. If only the peak discharge is needed, the method developed by the 
Natural Environmental Research Council of Great Britain (1975) can be used 
(see Appendix). 


38 












The best results can be obtained when the flood routing is done using the 
full equations of unsteady open channel flow (St. Venant equations) as shown 
in the studies by Huntington (1974) and Fread (1976). However, this type of 
analysis requires a much greater effort in terms of personnel time, is more 
expensive in terms of computer costs, and requires the hydrologist to have some 
background in numerical analysis techniques. 

Physical models have also been successfully used to evaluate the effects of 
encroachments. However, physical hydraulic models are costly to construct and 
test, and the time required to make a physical model study is usually long when 
compared to a computer modeling effort. When an existing model is available, 
studies of this type become competitive with numerical modeling studies and 
often produce more accurate results. 

4.2 Other Considerations 

As discussed above, in most cases when the flood plain encroachment causes 
no more than a 1 foot rise in the maximum water level, the corresponding increase 
in peak flow is usually minor. In addition, in practice the actual increase will 
usually be less because the full proposed encroachment often does not occur. 

As pointed out by Goddard (1979), the actual development on a flood plain is such 
that significant flood plain storage is still available even though the floodway 
fringe is developed. Also, the encroachment limits selected usually produce a 
water surface rise that is less than 1.0 ft. 

Even though the actual effects of encroachments on flood hydrographs will 
likely be less than the anticipated amounts, this cannot be relied upon in 
floodway analysis. However, it does provide a safety factor to help ensure 
that the actual water surface rise does not exceed the designed amount 
(say for example, the one foot used for FIA studies). 











4.3 Examples 

Some examples are presented to illustrate the circumstances for which peak- 
flow increases can be ignored and to indicate when additional analysis must be 
done. As discussed above, each case must be considered by itself, and specific 
numerical guidelines cannot be given. 

Example 1 - Steep River, Minor Encroachments 

A flood-plain study is being made for a river which has an average slope of 
five feet per mile. The encroachment limits which give a maximum rise in water 
level of one foot for the 100-year flood result in a removal of ten percent of 
the flood plain along a four mile reach of the river. The major floods in this 
river are snowmelt floods which last for several days. The question to be 
considered is: will there be significant changes in the flood peak in the 
areas downstream from this encroachment? 

In this case the removal of a relatively small volume from flood plain 
storage will have only a minor effect on the flood hydrograph, especially since 
the river slope is not flat (S = 0.001). In addition, the hydrograph is not 
sharply peaked, further reducing the attenuation effects. Making a detailed 
flood routing analysis is not warranted, therefore. 

Example 2 - Large River, Major Encroachments 

The Missouri River is joined by the Kansas River at Kansas City, Kansas 
and Missouri. Flood walls and levees along the top banks of the rivers 
are used to protect these cities against flooding. It was proposed to construct 
a series of agricultural levees along the river below Kansas City which would 
remove many square miles of flood plain area. Because major amounts of flood 
plain area would be no longer available, the effect of these levees on the 
flood conditions at Kansas City would be great. In the reach of the river that 
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was studied, the flood plain is from 1-1/2 to 11 miles wide and about 300 miles 
long. The levees limit the floodway width from a minimum of 3,000 feet to about 
5,000 feet. 

An evaluation of the effects of these encroachments was done using a physical 
model, the Mississippi Basin Model of the U.S. Amy Corps of Engineers (WES, 
1955). In this study it was found that for a design flow of over 500,000 cfs 
stages were increased by 3 to 5 feet along the river as a consequence of these 
levees. These stages were measured under steady flow conditions; the heights 
would have been somewhat less if they corresponded to maximum values resulting 
from a hydrograph routed through the model. 

In another test, a flood hydrograph with a peak discharge of 300,000 cfs was 
compared with a steady flow condition of 300,000 cfs. Stages for the steady 
flows varied from 0.5 feet to 1.5 feet higher than for the hydrograph flows. 

In general, simulations using steady flows will give higher flood stages than 
for unsteady flow conditions. 

Example 3 - Local Encroachment 

In a study of the effects of storage on flood magnitudes and stages for the 
North Branch of the Chicago River using the Hydrocomp Simulation Program (HSP 
model), the effects of a land fill 0.8 miles long by 0.4 miles wide were 
investigated (Hydrocomp, Inc., 1973). This is about 3 percent of the effective 
flood plain storage capacity. The channel slope is 3.8 feet per mile. This 
fill removes about 140 ac-ft of flood plain storage for the 100-year flood. 

At the land fill site, the computed peak discharge was increased 14 percent 


with an increase in the stage of 1.9 feet. At relatively short distances 
downstream the effects were negligible, however, with no apparent change in 
stage and less than 1 percent increase in peak discharge. 







Individual fill projects of this type usually do not displace large amounts 


of flood plain storage, and thus their influences are small, except in their 
immediate vicinity. The effects are cumulative, however, and several land 
fills of this type in combination may cause significant increases in discharge 
downstream. 
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APPENDIX 


FACTORS AFFECTING FLOOD WAVE ATTENUATION 

The description of flood routing procedures developed by the Natural 
Environment Research Council (1975) for use on British rivers provides an 
excellent source of information on factors affecting flood wave attenuation 
as well as general background on flood routing methods. The "convection- 
diffusion equation" is used for flood routing. The equation is written in 
a form which permits direct computation of flood wave attenuation along a river 
reach as a function of hydrograph and flood plain and river channel parameters. 
A brief summary of the assumptions used in its derivation, as well as 

limitations on its use, are discussed here. 

As with most flood routing procedures, the starting point is the St. Venant 
equations. The following assumptions can be made to permit simplification of 
the equations: 

1. River slopes are on the order of 0.001 or steeper. This implies that 
the momentum of the flow is governed primarily by the bottom and 
friction slopes. 

2. Acceleration and convection of momentum terms are small and can be 
ignored. 

3. Contribution of momentum from lateral inflows can be ignored. 

4. The volume of lateral inflow is small relative to channel flow. 

The resulting simplified equations of motion are: 



A-l 







where A is cross sectional flow area, t is time, Q is flow, x is direction 
of flow, q is lateral inflow per unit length of channel, S is channel slope, 
y is flow depth, n is Manning's roughness parameter, and R is the hydraulic 
radius. 

Further simplifications can be made if the following assumptions hold: 

1. The computation reach is short so that changes in flow, wave speed, 
and channel slope are small (changes are less than 10 percent, say) 

2. The flood wave speed is governed primarily by the main channel flow. 
The simplifications to the equations of motion enables routing of the 

peak flow without considering the full hydrograph. Such an equation is of 
interest here because it illustrates the factors which influence flood wave 
attenuation. 


At the flood peak aQ/at = 0. For a sufficiently short river reach this 
condition also holds along the entire reach. If lateral inflow is assumed to 
be negligible, it can be shown that the attenuation of the flood wave, Q* is 


given by: 



(3) 


where the subscript 1 denotes values associated with the peak flow at the 
upstream end of the reach. Thus, is the upstream peak flow, and is the 

downstream peak flow, so that Q* = a is an attenuation parameter for 

the reach, c^ is the mean flood wave speed for the reach: = l/T , where L 

is the reach length and T p is the time difference between peaks at the upstream 

2 2 

and downstream ends of the reach; d Q/dt is the second derivative of the flow 
with respect to time at the peak. This last term gives a measure of the 
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hydrograph curvature at the peak. The reader is referred to Natural Environ¬ 
ment Research Council (1975) for the details of the derivation and procedures 
for application to flood routing studies. This publication also contains some 
comparisons between this method and other routing procedures. 

The attenuation parameter a is a function of width and length of the flood 
plain plus river channel and of the river slope. Because the river width and 
slope are not constant, the routing reach can be divided into subreaches. 

If M subreaches are used 



where P is the plan area of the inundated flood plain and channel for the m-th 
m 

subreach, L m and S m are corresponding lengths and channel slopes. 

The curvature of the peak of the upstream discharge hydrograph is approximated 
by: 


».t ♦ «-» t • ”ii 

dt 2 (At ) 2 


(5) 


where Q 1 is flow at the peak, Q_ At and Q At are the flows at some time At before 
and after the peak respectively. 

The flood wave speed c (Q) is defined as the average speed along the reach 
at the peak discharge. In the derivation of Equation (3), c was based on a 
peak flow without attenuation. When there is attenuation, c depends not only 
on the observed flood wave speed L/T p , but also on the shape of the discharge 
hydrograph. It is recommended in the Natural Environment Research Council 
report that if the attenuation is greater than 10 percent (QVQ^ >0.1) 
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then c should be computed by the expression 

c = Ip- - b Q* (6) 

'p i 

The purpose of this discussion is to show what factors influence the 
attenuation of the flood wave. Under conditions which meet the assumptions 
listed above, the important factors as indicated by Equation (3) are: 

1. Hydrograph shape, especially curvature at the peak, 

2. Flood plain width, 

3. River slope, and 

4. Flood wave speed. 
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